Why Are Alkali Halide Solid Surfaces Not Wetted By Their Own Melt? 
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Alkali halide (100) crystal surfaces are anomalous, being very poorly wetted by their own melt at 
the triple point. We present extensive simulations for NaCl, followed by calculations of the solid- 
vapor, solid-liquid and liquid-vapor free energies showing that solid NaCl(lOO) is a non-melting 
surface, and that its full behavior can quantitatively be accounted for within a simple (BMHFT) 
model potential. The incomplete wetting is traced to the conspiracy of three factors, namely: surface 
anharmonicities stabilizing the solid surface; a large density jump causing bad liquid-solid adhesion; 
incipient NaCl molecular correlations destabilizing the liquid surface. The latter is pursued in 
detail, and it is shown that surface short range charge order acts to raise the surface tension because 
incipient NaCl molecular formation anomalously reduce the surface entropy of liquid NaCl much 
below that of solid NaCl(lOO). 

PACS numbers: 63.03.Cd; 68.08.Bc; 68.08.De 
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Molten alkali halides and their surfaces have long been 
studied experimentally [J Q and theoretically |^ y, Q IE 
0,0- Much less attention has been devoted to solid alkali 
halide surfaces at high temperatures, and especially to 
their wetting habit at the melting point. Yet, these crys- 
tal surfaces behave anomalously in that respect. Whereas 
most liquids would wet their own solid- vapor interfaces - 
the solid surfaces undergoing surface melting at the melt- 
ing point Tm (triple point wetting '^)- molten salts wet 
their own solid surface only incompletely. The external 
and internal contact angles 9 and of a partially wetting 
liquid droplet onto its own solid are connected to the free 
energies 7 of the solid- vapor (SV), solid-liquid (SL) and 
liquid- vapor interfaces by Young's equation [oj: 



ISV = ISL cos 4> + ILV cos ( 



(1) 



Incomplete, or partial, triple point wetting, with ~ 
but with 6 > Q and stable near T^, implies that some 
physical mechanism must be at work making 75^ < 
ISL + ILV- For liquid NaCl on NaCl(lOO) at = 
1074 K, a partial wetting angle ~ 48° is well estab- 
lished, with similar results holding for other alkali halides 
too. [ifll |ll| This angle is exceptionally large - and thus 
the corresponding self-wetting at the triple point excep- 
tionally poor - even when compared with strongly non- 
wetting metal surfaces, where e.g. 9 ~ 15° measured 
for liquid Pb/Pb(lll) 0|, or ~ 18° obtained by sim- 
ulation of liquid A1/A1(111) Should liquid surface 
layering Q be, as in the case of metals, ^2 be the 
culprit for the incomplete wetting of alkali halides too, 
the layering magnitude and its effects should be excep- 
tionally strong. However, all of the existing molten salt 
theory and simulations indicate the opposite, namely a 
soft, smooth, layering-free liquid-vapor interface 0,|Hl3- 
That leaves the partial wetting wide open for another ex- 
planation. In a more general context, it seems desirable 



to pursue a case study of solid-liquid- vapor coexistence 
at the triple point in a simple but realistic model that 
could be addressed at the fully microscopic level. We 
wish in particular to understand what may control wet- 
ting and adhesion in a specific and chemically unambigu- 
ous case such as that of a liquid with its own solid. For 
these reasons and in order to shed light on the under- 
lying physics, we undertook extensive simulations of the 
NaCl(lOO) solid surface, of the liquid NaCl surface, and 
of the solid- liquid interface around T,n. 

NaCl, our prototype alkali halide, was described by the 
classic Born-Mayer- Huggins-Fumi-Tosi (BMHFT) two 
body potential 'l6|. Polarization forces, though not neg- 
ligible 6], were sacrificed in the present context, where 
computational simplicity is essential to reach a unified 
description of all possible interfaces with very large sizes 
and simulation times. Eventually, as it turnes out, quan- 
titative accuracy in the description of NaCl interfaces 
seems anyway quite good in the BMHFT model. 

Bulk systems were first simulated by molecular dynam- 
ics (MD) at constant volume with cubic simulation cells 
comprising up to about 5000 NaCl units. Surfaces were 
studied with periodically repeated slabs - consisting of 
12-^24 planes with 64 NaCl units each - separated by 
100h~120 a of vacuum. We treated long range forces 
in full using a 3D Ewald summation. Despite the size 
and time limitations imposed by long range forces, great 
care was taken to run simulations long enough for a clear 
equilibration, typically 100-300 ps at Tm- In preliminary 
bulk simulations, we determined the equilibrium lattice 
spacing a{T) as a function of temperature at zero pres- 
sure. All subsequent simulations were then performed at 
fixed cell size, by enforcing a{T) at each temperature. 
The calculated room temperature aiTu) and linear ex- 
pansion coefficient, 5.683 A and 40.5 • 10~^ are 
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in very good agreement with experiment, 5.635 A and 
38.3 • 10"^ K-^ The bulk melting temperature - 
practically identical to the triple point temperature ow- 
ing to Clapeyron's equation - was established by sim- 
ulation of liquid-solid coexistence at zero pressure and 
found to be 1066±20 K (1074 K experimental). [13 This 
value is also in perfect agreement with 1064±14 K inde- 
pendently obtained for the BMHFT potential by Anwar 
et al [2d |. Our theoretical volume expansion and latent 
heat of melting are 27% and 0.2899 (experimental 26% 
and 0.2915 eV/molecule Q). Full details will be pre- 
sented in a forthcoming paper [l^ . 

In surface MD simulations, defect free NaCl(lOO) was 
first of all found to remain indefinitely solid and to- 
tally dry for all temperatures up to T™, and for even 
the longest simulation times ~ 1 ns. Moreover, in a 
metastable state, simulated NaCl(lOO) remained crys- 
talline even above Tm- By melting initially a few sur- 
face layers, and then observing in a subsequent canonical 
MD run that these liquid layers spontaneously recrystal- 
lized even above T^, |2j|, we established the existence of 
a multi-monolayer thick nucleation barrier against sur- 
face melting up to about 1115 K ~ Tm+ 50 K, then 
a thin monoatomic nucleation barrier until a 1215 K ~ 
Ts = Tm+ 150 K. Only at T^, well above r„, does sohd 
NaCl(lOO) become locally unstable and spontaneously 
melt. This metastable behavior was found to persist even 
in presence of surface defects, such as molecular vacan- 
cies [i^l or steps, at least up to 1115 K 0. Altogether, 
these results characterize NaCl(lOO) as a clear case of sur- 
face non-melting, 0, 0, 0| a prediction that deserves 
to be tested in experiment. For a short enough time, 
it should be possible for example to overheat NaCl(lOO) 
and other alkali halide surfaces substantially above Tm 
without melting them. 

We proceeded next to simulate the liquid NaCl surface, 
more precisely melted NaCl slabs. The (x,y) averaged 
liquid local density profile p{z) (Fig. ^) is confirmed 
to be remarkably smooth as found in previous stud- 
ies H,ISI3j devoid of layering unhke Al, Pb (Fig. n^,b), 
and much more diffuse (Fig. QJi) than even that of liq- 
uid Ar (Fig. n^). The nature of diffuseness of the NaCl 
surface is demonstrated by the simulation snapshot of 
Fig. showing very pronounced fluctuations [i^, in 
the instantaneous surface profile. This picture is sug- 
gestive of a low surface tension, high entropy surface, in 
apparent contradiction with the massive non-wetting of 
solid NaCl(lOO) by its melt, which is only favored by a 
sufhciently large jlv (Eq- O)- order to clarify the 
situation, we undertook direct calculations of interface 
free energies 7sy, Ilv^ and 75^. We obtained jlv as 
the surface stress of the simulated liquid slab (two inter- 
faces) through the Kirkwood-Buff formula: 
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FIG. 1: Simulated density profile of liquid surfaces 15| of (a) 
Pb, (b) Al, (c) LJ (Ar) in comparison with (d) p{Na) and 
piCl) profiles in NaCl at Tm\ solid line: average of the two. 
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FIG. 2: a) Simulation snapshot of the NaCl liquid surface at 
Tm- Notice the large thermal fluctuations, and some nearly 
molecular configurations highlighted in the outermost region; 
b) Coordination number N{Z) and density profile showing the 
smooth crossover from liquid {N = 4.6) to molecular vapor 
(TV = 1.3 dotted line). Density profile as in Fig. 0in units 
p/pbuik. 
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where: (a, /3) = {+, — ), ia and jP denote ions at site i or 
j, Z is the distance normal to the interface, L^, Ly are the 
{x, y) dimensions of the supercell and ay = -^{axx + cfyy) 
and (T_L — fJzz are the tangential and normal compo- 
nents of the stress tensor respectively. Here ( ) denotes 
a canonical average and X^i j O'^^^ all pairs of par- 
ticles. Moreover r^j = {xij,yij, Zij) is the interatomic 
distance, fapifij) is the force between atoms i and j. 
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FIG. 3: Calculated NaCl surface free energies, a) Solid-vapor 
7s v- Note the metastability up to about 150 K above Tm- 
Dots: effective harmonic approximation, b) Liquid-vapor 
-yLV- Experimental data from Ref. and references therein. 
'Ylv'- hypothetical liquid surface free energy calculated by set- 
ting A = for outer surface atoms with coordination below 
1.3 (highlighted in Fig. Hi). 



9ai3{j^ij', Z) are the Na-Cl, Na-Na, Cl-Cl pair correlation 
function measured in a slice centered at Z, Pa{Z) the av- 
erage density of ion a near Z and finally A is a parameter 
here equal to one, but inserted for later use. 

As shown in Fig. |Hb), 71, v so obtained compares well 
against experiment for example at T„j, ^lv — 104±8 
mJ/m^ against 116 mJ/m^. 

The free energy ^sv of solid NaCl(lOO) was obtained 
by thermodynamic integration jsv/T = 7sv,o/2o + 
j^^d{l/T')l\E{T'), where A£;(T) is the surface excess 
internal energy obtained by simulation, and Tq = 50 K 
was chosen as a convenient reference state (our simula- 
tions are classical and do not include quantum freezing) . 
The results of Fig. |3^a) show a relatively low "isvi char- 
acterized by a drop from 210 mJ/m^ at 50 K to ~ 100 
mJ/m^ at Tm- The drop confirms that the exceptional 
stability of solid NaCl(lOO) is significantly enhanced by 
anharmonicity above 600 K. To gauge the nature of an- 
harmonicity, we extracted from velocity- velocity correla- 
tions the T-dependent frequency spectrum of a slab and 
of a bulk sample with same number of molecular units. 
Treating both spectra as collections of harmonic oscilla- 
tors, two harmonic free energies can be calculated. Half 
their difference yields an effective harmonic 7gy"", whose 
milder drop with temperature, due to surface vibrational 
softening, is seen to recover about half the true solid sur- 
face entropy Ssv = —djsv /dT. The remaining half thus 
represents an additional stabilization of the solid surface 
by "hard" anharmonicity. 

Finally 7s l in the BMHFT model was calculated by in- 
serting in Eq. the calculated jsv, ILV and the contact 
angles 9 = 50° ±5°, ^ 0°, obtained by separate simu- 
lation of a liquid droplet deposited on NaCl(lOO). |3| 
We obtained in this manner 75^ — 36±6 mJ/m^. A 
charged hard sphere |2l] study along the lines proposed 
e.g. for neutral hard spheres niight be desirable in 
order to rationalize this result. Our relatively large 75L 
indicates poor adhesion of the liquid to the solid sub- 



strate, attributable to the unusually large difference of 
density and structure between liquid and solid. |l7l | 

While the result ^sv — ISL + iLV ~ 40 mJ/m^ now 
fully accounts for the non-melting of NaCl(lOO) and its 
incomplete triple point wetting, it does not yet clarify 
the role of the liquid NaCl surface in this context. That 
can be addressed in more depth by considering the calcu- 
lated temperature dependence of the surface tension 7^,1/ 
(Fig. . Strikingly, the temperature dependent drop of 
surface free energy shows a factor 2.6 lower surface en- 
tropy Slv = —d'jLv/dT of the liquid surface compared 
with that Ssv of the solid surface. This is contrary to 
naive expectations, based on the pictures of a relatively 
ordered, defect free solid surface and a very disordered, 
strongly fluctuating liquid surface. The presence of this 
"liquid surface entropy deficit" (SED) suggests some un- 
derlying short range order, that could also explain why 
the liquid surface tension, is here, surprisingly, as high 
as that of the solid surface. The surface order, if any, is 
clearly not layering: so what else could it be instead? 

The answer, as was foreshadowed earlier on 0], is that 
charge order, already important in bulk, plays an en- 
hanced role at the molecular liquid surface. If surface 
thermal fluctuations are indeed large, we find them re- 
vealingly correlated. For a Na"*" ion that instantaneously 
moves e.g., out of the surface, there is at least one ac- 
companying Cl~, also moving out; and vice versa. So 
while large surface fluctuations smear the overall liquid 
vapor density profile, bridging gently between the liquid 
and the vapor, (Fig. ^i) the two-body correlations, de- 
scribed e.g., by the the Na-Cl pair correlation function 

g^i (r), or by its integral, the ion coordination number 

N , drop from values typical of the bulk liquid at Tm, to 
the nonzero value of the molecular vapor, instead of zero 
as in the LJ liquid. For a more quantitative understand- 
ing, we calculated a locally defined charge coordination 
number: 

N±{Z) = r^'^ \dZ'p{Z') [ d^r g+_(r;Z') 

(3) 

where = 4 A corresponds to the first local minimum 

of (r), and ^ is a small interval. Starting with the 

ideal Ns = 6 of the solid, we have Nl — 4.6 in the bulk 
liquid at T,„. Moving across the liquid- vapor interface 
N{Z) drops continuously from 4.6 downward (Fig. [SJs). 
Even if simulation statistics is lost in the vapor, N{Z) 
must remain and remains for large Z clearly larger than 
or equal to Ny — 1.3, the value appropriate for NaCl 
vapor (which at Tm consists for 69% of NaCl molecules, 
and 31% of dimers H^). The larger the coordination 
number of atoms in the interface region, the less their 
configurational entropy, the higher the surface tension. 
Hence incipient molecular order could be the reason for 
SED of liquid NaCl. 

For a test of this idea, we repeated the surface tension 
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calculation of liquid NaCl by only slightly and artificially 

altering in Eq. (0) the value of correlations or, which 

amounts to the same, of the forces acting among Na and 
CI for the (extremely small) fraction of outermost surface 
atoms whose coordination number iV ^ 1.3. Since the dy- 
namics, and thus the internal energy remain untouched in 
this way, the contribution of these configurations to the 
surface stress gives a direct measure of the negative sur- 
face entropy contribution to the surface tension caused 
by molecular short range order at the liquid-vapor inter- 
face. We first identified in simulations the surface Na 
and CI ions by means of a simple algorithm. We then ex- 
tracted their mean electrostatic potential value 1^, and 
established a roughly linear connection between atom po- 
tential and coordination E^ . Finally, we modified 
the Coulomb part of the charge contribution in Eq. Q 
of surface Na"*" ions in the form A — Q{Vq — Vi) where & 
is the step function and Vb = — 6.99 eV is the value that 
cancels correlations for Na+ ions with N < 1.3. Though 
representing an exceedingly small fraction of the surface 
atoms (Fig.|3i), removal of the surface stress contribution 
by these molecularly paired Na and CI ions yields a con- 
siderable surface tension drop from — 104 mJ/m^ 
to jIy = 53 mJ/m^ (Fig- The increased temper- 
ature slope \dj1y/dT\ exactly matches the calculated 
drop from jlv to 72y, confirming that it corresponds 
to a purely entropic gain - the removal of some of the 
SED through cancellation of molecular surface correla- 
tions. Since now + jsl < isv, (the equivalent cor- 
rection to and 75^ is utterly negligible) recovery of 
this surface entropy actually suffices to provoke complete 
instead of partial wetting of NaCl(lOO) at the melting 
point. 

In conclusion we found that BMHFT potential simu- 
lations and the related surface thermodynamics explain 
quantitatively the incomplete wetting of NaCl(lOO) by 
liquid NaCl at the triple point. Three elements, namely 
the exceptional anharmonic stability of the solid (100) 
surface, the poor adhesion of the liquid onto the solid, 
and a liquid surface entropy deficit caused by incipient 
molecular short range charge order, all conspire to give 
rise to this phenomenon in the BMHFT model of NaCl 
- and most likely also in real alkali halide surfaces. 

Experimentally, it should be possible to demonstrate 
the overheating of NaCl(lOO) and other alkali halide sur- 
faces, perhaps in the same way as in metals j27j (though 
here the high vapor pressure at suggests using tech- 
niques not relying on ultra-high vacuum). The poor ad- 
hesion of the molten salt onto its own solid should be 
detectable in nucleation. The short-range correlations 
described at the surface of molten salts could possibly 
become accessible spectroscopically. Also, if these sur- 
face correlations could for example be altered by external 
means, e.g., by electric fields, the wetting angle should 
change accordingly. Finally the general possibility that 
some form of short range order at the liquid-vapor in- 



terface might affect the surface tension by reducing sur- 
face entropy could prove of wider relevance to other com- 
pound and molecular liquid surfaces. 
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